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Abstract 

We develop a theory based on Bohmian mechanics in which par- 
ticle world lines can begin and end. Such a theory provides a realist 
description of creation and annihilation events and thus a further step 
towards a "beable-based" formulation of quantum field theory, as op- 
posed to the usual "observable-based" formulation which is plagued by 
the conceptual difficulties — like the measurement problem — of quan- 
tum mechanics. 

PACS numbers 

03.65.Ta (foundations of quantum mechanics) 

42.50.Ct (quantum description of interaction of light and matter) 

1 Introduction 

A formulation of nonrelativistic quantum mechanics based on objectively ex- 
isting particle positions and particle trajectories, now usually called Bohmian 
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mechanics, was proposed fifty years ago by David Bohm (and even earlier 
by Louis de Broglie); see j3] for a recent overview. Today there remain two 
big challenges for this approach: to form a relativistic version, and a version 
suitable for quantum field theory (QFT). Here we shall address the latter. 
We describe a general and particularly natural way of extending Bohmian 
mechanics to QFTs, explicitly giving the equations for a simple example. 

Bohm himself proposed [21 p. 230] that Bohmian mechanics should be 
extended to QFT by means of the incorporation of the actual field configu- 
ration, guided by a wave functional (the state vector). In contrast, John Bell 
proposed a model [TJ p. 173] in which, instead of the field configuration, the 
local beables are the fermion numbers at each site of a lattice discretizing 3- 
space. We argue that it is instructive to modify Bell's proposal in two ways, 
and thus get a similar but even simpler theory with a direct connection to 
Bohmian mechanics. Below, we give an explicit example of such a theory for 
a particularly simple Hamiltonian. 

Bell's model contains no beables representing the bosonic degrees of free- 
dom (such as radiation), neither an actual field, nor actual particles, nor 
anything else; the existence of a radiation part of the state vector is relevant 
only to the behavior of the fermions. This is certainly consistent and empir- 
ically irrefutable, but it is neither necessary (as our example below shows), 
nor even a natural view. 1 In the model proposed below, the bosons have the 
same status as the fermions: they are particles, described by their positions. 
This, of course, should not be regarded as discouraging consideration of the 
approach based on actual field configurations. 

The other deviation from Bell's proposal is the replacement of the lattice 
by continuous space. The lattice was introduced in the first place for the 
purpose of providing an effective ultraviolet cut-off and thus a well-defined 
Hamiltonian. Such a cut-off, however, can also be realized by smearing out 
the interaction Hamiltonian through convolution with, say, a sharply peaked 
but bounded function ip. The continuum analogue of the particle number 
at every lattice site is the position of all particles in ordinary space, with 
the total number of particles N(t) possibly varying with time. However, in 

1 After all, if the role of the wave function is to guide particles, it would seem that 
there should be as many particles as there are variables in the wave function. In addition, 
the Hilbert space of the radiation degrees of freedom (unlike the one for quark color) is 
not a purely abstract Hilbert space, but is related to space-time points (via creation and 
annihilation operators, for example), a fact that would seem surprising if the state vector 
were not related to space-time objects such as particles, strings, or fields. 
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the model we propose the particles follow Bohmian trajectories, except when 
there is particle creation or annihilation. (For a discussion of how Bohmian 
mechanics arises from a lattice model (in the absence of interaction) in the 
limit of vanishing lattice width, see [HIE]-) 

Our proposal profits, we believe, from making this contact with Bohmian 
mechanics, since then every argument for taking the Bohmian trajectories 
seriously also provides some support for the proposal. Bohmian mechanics 
also profits from this contact because the Bohmian trajectories can then be 
taken seriously even in the framework of QFT. 

Our moving configuration Q(t) is constructed in such a way that it is 
random with distribution at time t equal to p(t) = \^(t)\ 2 , where \l/(t) is the 
(position) Fock space representation of the state vector ^(t)) at time t. (For 
a model in which there are only bosons, the n-particle component of ^ is 

¥ n \ Xl , ...,x n ) = -^= (0\a( Xl ) ■ ■ ■ a(x n )\V) (1) 
Vn! 

where |0) is the Fock vacuum and a(x) = (27r) _3//2 Jd 3 k e lk x is the bo- 
son annihilation operator at position a?.) In particular, the probability that 
N(t) = n, i.e., that there are n particles at time t, equals the integral of p(t) 
over K^ 1 ' = (M 3 ) ra (suitably symmetrized), the n-particle sector of configu- 
ration space — i.e., the L 2 norm of the projection of onto the n-particle 
subspace of Fock space. A nontrivial superposition of quantum states with 
different particle numbers leads to a probability distribution over different 
particle numbers, among which only one is, of course, actually realized. 

2 Configuration Jumps 

The Hamiltonian for QFTs is typically a sum of terms, each of which yields a 
contribution to the motion we wish to propose. It is quite generally the case 
that H = Hq + H irit where the free Hamiltonian H corresponds naturally 
to a deterministic motion — in the model considered here that of Bohmian 
mechanics, and in a relativistic model for example the Bohm-Dirac motion 
[21 p. 274] — given by a (time-dependent) velocity field v = v^^(q,t), on 
configuration space K = U n K^ n \ that defines the "deterministic part" of 
the process, while the interaction term H mt corresponds to random jumps, 
a continuum version of the process proposed by Bell |1, p. 173], defined by 
jump rates a = cr(q', q, t) = er*(*)(g', q) for a transition from q to q' at time t. 
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[This means that when the actual configuration Q at time t is q, then with 
probability (density, with respect to q') a(q', q, t) dt, Q will jump from (very 
near) q to q' in the time-interval (t, t + dt).} 2 

The relevant continuum analogue of Bell's jump rates for H = Hi nt , 
equations (6-8) of [H p. 173], is 



where for particles with spin the two products in the numerator and the 
product in the denominator are local spinor inner products, and where we 
have used the notation A + = max (A, 0) for the positive part of A 6 R. This 
will typically be well-defined (as jump rates), since the kernel (g|if int |g') of 
the (cut-off) interaction Hamiltonian of a QFT should involve nothing worse 
than ^-function singularities. 

The complete process, corresponding to the total Hamiltonian, is then 
given by the deterministic motion with velocity v, randomly interrupted by 
jumps, with rate a, after each of which the deterministic motion is resumed 
until it is again interrupted. As a function of t, each realization Q(t) is thus 
piecewise smooth. At the end of a smooth piece, Q jumps to the starting 
point of the next smooth piece. What is stochastic about Q are the times 
at which the jumps take place, and the destinations of the jumps. The 
probabilities for times and destinations are governed by the wave function. 
Q{t) is a Markov process. 

As described, our process is so designed as to directly imply the following 
equivariance theorem: IfQ(to) is chosen at random with distribution |\I/(to)| 2 , 
then at every later time t > t Q , Q(t) is distributed with density |\l/(t)| 2 . (This 
can also be explicitly checked by comparing the equation for d(^^)/dt as 
implied by the evolution equation of the quantum state, ihd^/dt = H^fr, 
and the master equation for the distribution of Q(t), which reads 

—p{q, t) = -V- (p{q, t) v{q, t)) + J^dq'(p{q', t) a(q, q', t)-p{q, t) a(q f , q, t)) . ) 

The choice of jump rates that will make the |\1/| 2 distribution equivariant 
is not unique. But the 'minimal' jump rates are unique, and these are the ones 

2 One can regard the jumps as having two ingredients: A jump occurs with total rate 
er(q, t) — J a(q' , q, t) dq'; when a jump does occur, the destination is randomly chosen with 
distribution a(q' , q, t)/a(q, t). 
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we have chosen. Details on the choice of jump rates as well as an elucidation 
of the general mathematical structure underlying equivariant processes will 
be presented in [3]. 

Corresponding to the standard interaction terms in QFT, the possible 
jumps are very restricted: there are only changes in the particle number by 
±1, and possible types are (i) appearance or (ii) disappearance of a particle 
(while the others remain at their positions) or (hi) replacement of a particle 
by two others or (iv) the reverse of this. Type (i) is appropriate for, e.g., 
photon emission, (ii) for absorption, (iii) for pair creation and (iv) for pair 
annihilation. (It follows that there is no discontinuity in the individual par- 
ticle world lines, in spite of the discontinuity in Q.) In our explicit model, 
only types (i) and (ii) occur. 

3 An Explicit Model 

We now present an explicit theory. It is based on a "baby" QFT taken from 
[HI p. 339] and [5], containing two species of particles, which we simply call 
electrons and photons. Electrons, whose number stays constant, emit and 
absorb photons. To keep things simple, we employ nonrelativistic dispersion 
relations (as everything here is nonrelativistic) for both electrons and pho- 
tons, and give the photon a positive (rest) mass. In addition, we ignore spin 
and polarization, and, of course, smoothen the interaction Hamiltonian. 

In held theoretic language, we have a bosonic held <p(x) = a\x) + a{x), 
with a) and a the photon creation and annihilation operators [cf. equations 

and (j5J)] and a fermionic held ifj(x), and the Hamiltonian is the sum 

H = H F + H B + H mt 

= (l^/^V^V^) + m mB) j^Va^)Va (x) 

+g Jd 3 x^(x)4> v (x)ilj(x) 

where 4>ip(x) = J d?y(ip(x — y)a){y) + Tp(x — y)a(y)) is the cutoff bosonic 
held, m-p and me denote the mass of the electrons and photons, and g is 
a real coupling constant. H commutes with the fermion number operator 
N F = $d z x^(x)ijj(x). 

Since the fermion number is conserved, we give it a fixed value N. The 
configuration space (changing notation slightly from before) is K = Um=o 
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where m is the photon number and 

K (m) ._ ( R 3)7V x ( R 3)m ? ^(0) = ( R 3)7V _ 

We will denote the electron coordinates by x := := (x 1 , . . . ,x N ) and 
the photon coordinates by y :— = (y 1 , . . . , y m ). The full configuration 
is thus given by q = (x,y) and, more explicitly, by = (x,y^). 

In terms of the wave function \l> [the (position) Fock representation of the 
quantum state], the Hilbert space inner product is 

oo „ 

($|M>> := / d 3N xd 3m y $(x,y)*(x,y), 
and the contributions to the Hamiltonian assume the form 



1m\ 
h 2 



2m B 3 
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and 



N N 



tfmt = g J2 M*) = 9 E + M^)) (3) 

i=i i=i 

with creation and annihilation operators and acting on Fock space in 
a smeared-out form: 

ajp(x) = j ' d 3 uTp(u — x) a(u), 
(at( W )*)( g M) = _L£a(yi -„)*(?■), ( 4 ) 



where 



(a(u)*)(g< ro )) = Vm+7#(g (m) ,i*). (5) 

aj,(aj) creates a new photon in state if centered at x e ffi 3 (which will be the 
position of an electron), and a^(x) annihilates a photon with "form factor" 
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<p at x. 3 Here (q^ m \ u) is the configuration with a photon at u added to q^ m ' 
and g- 3 is the configuration with the j-th photon deleted from q( m \ Thus 

N m 

Jm)\ _9_ 



i=l 3=1 

+ gV^+TYl /^W-* < )*(? M ,l/ , ) • 

1=1 ^ 

\1/ satisfies the Pauli principle, i.e., it is symmetric in the photon variables 
and antisymmetric in the electron variables. 

We now turn to the particles, described by the actual electron configura- 
tion X and the actual photon configuration Y. The deterministic part of the 
motion Q{t) = (X(t),Y(t)), corresponding to H F + H B , is the usual Bohm 
motion 

x = v ^ Q) ■= ^ Im m ' (7) 



vi » ^ ^ t (^W)(g) /on 
F = := ^ Im g(Q) (8) 



and it follows from (J2J) and (JHj) that only two kinds of jumps occur, with rates 
as follows: 

• The j-th photon vanishes, while all the other particles stay at their 
positions. Thus q' = qi and the jump rate from q = q^ to q' is 

a{q > q) = V^[-^ Im 9® J • (9) 

• A new photon appears at y', while all the other particles stay at their 
positions. Thus q' = (q^ m ',y') and the jump rate from q = q^ to q' 
(a density in the variable y') is 

<r{q,q) = —Vm+l[-}^lm — j . (10) 



3 If the form factor <p is square-integrable (as we assume) it provides an ultraviolet 
cutoff; it can be regarded as determining the effective range of a electron's power to create 
or absorb a photon. 
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We reiterate that these choices guarantee that the process obeys the equiv- 
ariance theorem mentioned in the previous section. 

Note that, contrary to what might have been expected, the creation op- 
erator in the interaction Hamiltonian (J3j) corresponds, according to (J2J), to 
the annihilation rate, while the annihilation operator corresponds to the cre- 
ation rate. There is less in this than meets the eye: The correspondence is 
an artifact of the way (0) is written, and if (J2J) had been written equivalently 
as 



the correspondence would have been reversed. Note also that if is supported 
by a 5-ball around the origin of M 3 , then photons can be created or annihilated 
only in the (^-neighborhood of an electron. 4 

Our model is translation invariant (unless one introduces an additional 
external potential), time translation invariant, rotation invariant provided 
if is spherically symmetric, time reversal invariant provided <p is real-valued, 
and gauge invariant provided one introduces an external vector potential — in 
all derivatives acting on electron coordinates — and an external scalar poten- 
tial into Hp. Galilean boost invariance fails, but not because of the Bohmian 
variables; owing to the interaction term (jB|). it fails, in fact, for the evolu- 
tion of the quantum state. The reason is that, roughly speaking, a photon 
gets created with wave function ip which cannot be Galilean invariant. Al- 
ternatively, one could say that under boosts with velocity u the form factor 
ip(y — x) must be replaced by exp (im B w • y/h) ip(y — x). Note that even 
with the cutoff removed, i.e., for ip(y) = S(y), the quantum dynamics is not 
Galilean invariant. 

If we add to if a suitable confining potential (and (p is real-valued, and 
the rest energy of the photon is made positive), it possesses a unique ground 
state [7j, and this ground state (like every nondegenerate eigenstate) is real 
up to an overall phase. Thus in this state the jump rates ©, (fTUJl as well 
as the velocities (|7Jl. (JHJ) vanish identically — nothing moves. Surprisingly, 
perhaps, the Bohmian particles do not perform any "vacuum fluctuations" . 

4 One might suspect that if S were small, there would be only a small probability 
(perhaps of order <5 3 ) that a photon will come closer than S to an electron. But equivariance 
implies that this is not so: if a certain amount of |\I>( m )| 2 flows to the sector i^( m_1 \ then 
the probability for a photon to be annihilated is just as large. 
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As in Bohmian mechanics, disentangled subsystems are governed by the 
same laws as for the whole. More precisely, if $ and $ are two wave functions 
from Fock space having disjoint supports S$ and in physical space M 3 , 
$ y& defines another Fock state $ ^ after suitable symmetrization; if the 
distance between S$ and is at least the diameter 5 of the support of <p, 
then in a system with wave function $ ^ the particles in S<& and those in 

will not influence each other, each set moving independently according 
to the corresponding 'factors'. This is a consequence of the fact that the 
velocities and jump rates are homogeneous of degree in the wave function. 

While the other particles keep their positions when a photon is created 
or annihilated, their velocities may change discontinuously because (J7|) eval- 
uated at the destination may differ from (J7J) evaluated at the point of de- 
parture. As a result, the world lines of all (possibly distant) particles, if 
the particles are entangled, will have kinks at the times of particle creation 
or annihilation. These kinks will however not be visible in, say, a cloud 
chamber since the necessary entanglement is destroyed by the decoherence 
of the tracked particle with its environment, caused, say, by the particle's 
interaction with the vapor. 

4 Removing the Cutoff 

Removing the cutoff is of course problematical, which is why the cutoff was 
introduced in the first place. However the problems arise primarily from the 
evolution equation of the wave function, not those of the Bohmian configura- 
tion, equations (J7J) through (fTUj). Suppose that a family of <^'s is parametrized 
by A G 1, and that as A — > oo, (p A (y) approaches S(y). Then the limit 
A — > oo corresponds to "removing the cutoff". Unfortunately, the corre- 
sponding H will not converge in any reasonable sense. But there exist 
numbers E A tending to infinity [3] such that H A — E A does converge in a 
suitable sense, and the evolution of the wave function is well-defined. This 
seems completely acceptable. We do not know whether the corresponding 
process Q A approaches (given a fixed initial wave function) a limiting process 
Q°°. To decide this requires a careful mathematical study, but at least we 
see nothing precluding this possibility: the velocity law is not affected by the 
cutoff, while absorption might become deterministic in the limit A — > oo and 
occur whenever (and only when) a photon hits an electron. 

Other Hamiltonians, more sophisticated ones, are more problematical, 
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and do not possess a limit as A — > oo, even after subtracting an "infinite 
energy" . In some cases this is due to the creation of a large (average) number 
m A of photons that goes to infinity with A. On the other hand, it is not clear 
that removing the cutoff is desirable or relevant. That is, there might exist 
an effective UV cutoff in nature, just as there is an effective IR cutoff (the 
finite radius of the universe). 

Be that as it may, if the unitary evolution on some Hilbert space does 
not survive the limit A — > oo, we face a problem, one that seems particularly 
bad for Bohmian theories, which so heavily rely on the wave function and, 
consequently, its having a well-defined unitary evolution. But appearances 
are misleading here. From a Bohmian viewpoint, the basic variable, bearing 
all the physical implications of the theory, is the configuration Q, whereas 
\1/ and H are only theoretical objects whose purpose is to generate a law of 
motion for Q. And, indeed, a law for Q might arise as a limit A — > oo of the 
law induced by \I/ A and if A ; it might be the case that while \1/ A and H A do not 
have a limit, the time evolution of Q is well defined in the limit. After all, this 
is precisely what occurs when one consider the limit H — > of nonrelativistic 
quantum mechanics: while the Hamiltonian and wave function do not have 
any sensible limit, the law for Q does! 

5 Determinism 

We close with a remark on (the lack of) determinism. It may seem surprising 
that we abandon determinism. Was not the main point of hidden variables 
to restore it? Actually, no. What was important was to provide a clear 
and coherent account of quantum mechanics. The simplest such account, we 
believe, is provided by Bohmian mechanics, which happens to be determin- 
istic. And the simplest such account of QFT seems to be of the sort we have 
presented here, which is stochastic. 
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